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Abstract 

Early universe equations of state including realistic interactions between constituents are built 
up. Under certain hypothesis, these equations are able to generate an inflationary regime prior 
to the nucleosynthesis period. The resulting accelerated expansion is intense enough to solve the 
flatness and horizon problems. In the cases of curvature parameter k equal to or +1, the model 
is able to avoid the initial singularity and offers a natural explanation for why the universe is in 
expansion. All the results are valid only for a matter-antimatter symmetric universe. 
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I. INTRODUCTION 



The gravitational field, as described by General Relativity, couples to all types of energy: 
rest masses, kinetic terms and interaction terms. Relativistic cosmology is, in consequence, 
deeply concerned with such sources. The kinetic terms and the rest-masses are currently 
taken into account by assuming ideal cosmic fluids constituted by ultrarelativistic and/or 



non-relativistic matter. Solutions of this type are found in standard texts 



2| and - when 



multi-component fluids are considered - in several papers, e.g., Interaction terms are 
commonly used only in perturbative models. In fact, using the Boltzmann equation in a 
Friedmann-Robertson- Walker (FRW) background, the inhomogeneities both in the cosmic 
microwave background (CMB^ 



with the observational data 



and in the matter content {4, 5| are studied, for comparison 
^, l7|. Notwithstanding, interactions are not considered as direct 
sources of gravitation in this line of research j^. 

Of course, some well-known proposals do consider interaction processes as direct sources 
of gravitation. They are usually related to the accelerated expansion regimes: present- 
day dynamics or inflation Nevertheless, these theories do not actually consider the 
fundamental interactions (electromagnetic, weak and strong) between particles in the source 



constituents. Indeed, in the standard inflationary approaches 



11 



an accelerated 



expansion is obtained through self-interaction processes of scalar inflaton fields (pi (x). This 
self- interaction is chosen so as to produce just those features which are necessary to describe 
the early accelerated regime [lOl. The phenomenological explanations for the present-day 

nnnnn 

acceleration include, among others, (i) the quintessence models [U, ll5|, ll6|, ll7|, ll8[ , which 
roughly follow the same lines of the inflationary theory; (ii) models of matter and dark- 
energy unification via Chapligyn-like equations of state (EOS) [l^ or through equations of 



20 



211 ]: (in) mode 



the Van der Walls type 
neutrinos to a quintessence scalar field 
energy conservation equation {23! . 



24|. 



s of mass- varying- neutrino type, which couple 
22|; (iv) models introducing interactions in the 



A formal procedure has been recently proposed for including the (fundamental) interac 



tions as direct sources of gravitation in the cosmological context 



. Imported from equilib- 



rium statistical mechanics, t 



nis 



state for both relativistic 



25 



brmalism allows the construction of realistic equations of 



26| and non-relativistic systems 



27 



28 



Our objective here 



is to find primeval cosmic fluid EOS taking into account physically realistic interaction pro- 
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cesses between the constituent particles, in addition to their kinetic and rest-mass terms. In 
particular, we shall examine their effect on the scale factor evolution. It will be shown that 
under certain hypothesis and approximations an early accelerated regime can be obtained 
as a consequence of these interacting processes. 

The idea of building realistic equations of state considering interaction between ele- 
mentary particles in the primeval universe is not new. Actually, during the late 1970's 
and the be ginn ing os the 80's a series of papers by Bugrii, Trushevsky and Beletsky 
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30 



31 



32l | discussed the construction of the high-energy EOS and their application 



to the pre-nucleosynthesis universe.^ Nevertheless, their treatment is different from the one 
developed here as we will make clear some pages ahead. 

The paper is organized as follows. Section [TTl reviews the general features of the early 
universe in its standard presentation , [lO| . Section IIIII presents some results from 



equilibrium statistical theory of interacting systems. Specifically, the coefficients appearing 
in the perturbative fugacity expansions are expressed in terms of the scattering matrix 
operator 5*. In addition, the matrix 5*2 describing the two-particle scattering is associated to 
the experimentally observable phase-shifts. The goal of Section IIVI is to construct realistic 
EOS for the pre-nucleosynthesis universe and discuss the hypothesis and approximations 
undertaken. In Section |Vl the more direct cosmological consequences coming from these 
equations are examined, including the effect of driving an accelerated expansion (inflationary 
era). Section IVll contains some final comments. Details of a too technical nature, as well as 
some data, have been relegated to appendices. 



II. GENERAL FEATURES OF THE PRE-NUCLEOSYNTHESIS UNIVERSE 

The period of the Universe evolution going under the name "Early Universe" covers many 
different and physically significant events. Indeed, it is usual to consider both the matter- 
radiation decoupling [kT^ ~ 0.1 eV) and the electroweak transition [kT^ ~ 100 GeV) as 
belonging to that period. This work is concerned with a specific interval within it, namely 
the pre-nucleosynthesis period (PNS). We shall define PNS as the period immediately before 
the nucleosynthesis of the light elements, when the nucleons (protons and neutrons) are in 

^ The authors are thankful to an unknown referee for caUing their attention to these works. 
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thermodynamical equilibrium with the rest of the cosmic fluid, i.e., kT^ > 20 MeV. So, 
in principle, any energy value larger than 20 MeV belongs to PNS; nevertheless, it will be 
enough to restrict our working frame to energies of a few hundreds of MeV. In consequence, 
PNS will in what follows actually mean the interval 20 MeV < kT < 300 MeV. 

Present-day universe is dominated by a dark energy component A plus a non-relativistic 
contribution formed by ordinary baryons b and dark-matter dm; the ultrarelativistic compo- 
nents - photons 7 and neutrinos u - are quite negligible. As one goes back in time, however, 
this situation changes drastically. Indeed, the evolution equations for ultrarelativistic and 
the nonrelativistic matter [l, 2! in terms of the expansion parameter, 

P-f,^ — P-fO,uO y~ J ^^'^ Pb,DM — PbO,DMO y~ J ' 

together with the Friedmann equations, 

SttG k a f d\ ArcG , ^ A 

and - = ^(p + 3p) + -, (2) 



Sa^S V*^/ 3 3 

show that the ultrarelativistic contribution overcomes that of the other components at en- 
ergies corresponding to kT^ ~ 1 eV .'^ The cosmological period dominated by ideal ultrarel- 
ativistic particles is usually called the radiation era. 

To find out the energetically relevant constituents during the PNS period we must notice 
that, as we turn to the past and the energy kT increases, pair-production processes become 
more and more frequent, generating a large variety of particles. The main scenario is that 
of ultrarelativistic particles generating non-relativistic ones. If we restrict ourselves to a few 
hundreds of MeV , the relevant particles are: 

• fundamental bosons: photons (7); 

• leptons: electrons (e~), positrons (e"*"), muons antimuons (/i'*'), electronic and 
muonic neutrinos (z/g, z/^) and electronic and muonic antineutrinos {ue,i^fi)', 

• hadrons: pions (tt^, tt", 7r°), kaons {K~^, K~ , K^, K^), nucleons (p, n) and antinucleons 
(p, n). 
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We will see in Section HVl that this statement must be qualified when the interaction processes are taken 
into account. 
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Concerning this list, it is important to emphasize that: (i) the particles taken into account 
are those with rest-mass bellow 1 GeV; and, (ii) the hadrons considered are only those which 
are stable by the strong interaction. Furthermore, we assume matter-anti-matter symmetry. 

Among the particles cited above, those that are ultrarelativistic at 20 MeV are 7, e^, 
Ue, z/^, z/r and z/^; the remaining ones are non-relativistic. By arguments given in classical 



texts 



5|, the chemical potentials involved in the PNS reactions are zero, i.e., /ij ^ for 
every species (i = 7, e^, /i^, ti"^, K^, ...). 

Thermal equilibrium is warranted as long as every reaction rate Ti{kT) of each given 
component i with all the other particles is much larger than the universe expansion rate. 
This last, in turn, is measured by the Hubble function H{kT). Therefore, when the condition 

r(T) > H{T) (3) 

is satisfied, thermodynamical quasi-static expansion holds. In fact, this may be verified for 
each variety by remembering that T = na \v\, where n is the target-particle density and a \v\ 
is the interaction cross-section times the relative velocity of the particles. The presence of 
thermal and chemical equilibria in PNS is extremely convenient, as it enables us to adopt, 
in the comoving reference frame, the usual statistical mechanics on the E'^ manifold (the 
phase space measure is given simply by j^^^n^^^^^P)- 
The energy density turning up in ([2]) is l|, [Sj 

9 Qtt^ 

Pur = P7 + Pe- + Pe+ + Pu. + Pv^ + PPe + = 2 = '30'*'^"^''^ 
and can be used to obtain the Hubble function as a function of kT: 



H,,T,J^.J'-^^^, (5) 

\a J V 5 mpi 

where nipi = = 1.221 x 10^^ MeV . The relation between the scale factor and the 

energy kT is determined by comparing to ([1]): 

kT = ^. (6) 

A is a constant which depends on the present-day values To and Oq. Substituting this 
last equation into ([5]) and solving the resulting differential equation - imposing the initial 
condition a{t = 0) = - leads to 



From ([6]) follows then the relation between the energy kT and the cosmological time t: 



The last four equations determine the evolution of the primeval universe in a radiation- 
dominated era. 



III. STATISTICS OF INTERACTING SYSTEMS 



It was discussed elsewhere |8[ how to construct equations of state for an interacting system 
of particles within the standard ensemble formalism of statistical mechanics, and how this 
could be applied to a simple model relevant to cosmology. In what follows, we will present 



a brief review of the information needed here. We also indicate Refs. 
for further information. 

The grand canonical partition function is expressed as 



25|, |28j, y and |35| 



eiz,v,T) = J2QN{v,T)^ 



.N 



(9) 



N=0 



where V is the volume, T is the temperature, z = e'^/^^ is the fugacity and Qn{V,T) is the 
N-particle canonical partition function. The chemical potential fi = fi^^ + m is composed 
by the rest mass m plus the non-relativistic chemical potential fi'^^. This last is that usually 



found in standard statistical mechanics texts, e.g. [27|]. Here, we are concerned only with 
one-component systems. For the cases of more constituents, see Appendix |X1 

The thermodynamical quantities - pressure p, energy density p, numerical density n, etc. 
- are related to the grand canonical potential 



1 °° 
n{z,T)^-\neiz,V,T) = Y,bNz'', 



(10) 



N=l 



written in terms of the cluster integrals h^- The thermodynamical limit has already been 
taken in Eq. (ITUi) . so that the h^q are functions of the temperature solely (the V dependence 
disappears). According to Dashen, Ma and Bernstein 25|, the cluster integrals are calculated 
from the ^-matrix as follows: 



N y 



-Tr I AS-^^S 
Ani \ dE 



dE, 



(11) 



CJV 



where j3 = 1/kT, gj~^ counts the degeneracy coming from internal degrees of freedom, A 
symmetrizes (antisymmetrizes) the bosonic (fermionic) states, 5* is the scattering matrix 
operator and fej^^ is the cluster integral of the non-interacting quantum system. The subscript 
Cat represents all the A^-particle connected diagrams'^ in which the interaction occurs at least 
once. In ( TTTl) . it was used the short-cut 

By using (ITO!) p, n and p are determined as series in the fugacity: 

( T'^ °0 oo rj, 

^ = J:6..-; n[z,T) = Y^Nh,z-- p(., T) = - J] (13) 

N=l N=l N=l ^ 

This is the parametric form of the equations of state. 

An alternative description is given by the pressure p{n, T) and the energy density p(n, T) 
written in terms of n (z, T). They are obtained by inversion of series n (z, T) and substitution 
into p{z,T) and p{z,T). This results in the virial expansion: 

'-^-pM ^J^^p,iT)n', (14) 

with ai and q representing the virial coefficients for the pressure and the energy density. 
These coefficients are completely determined by the b^. Appendix El presents the explicit 
forms of ai and q for a two-component system. 

In practice, it is not possible to really sum any of the series fl 13111 41) . Therefore, the choice 
of using p and p in terms of 2; or n depends on the perturbative characteristics of each series 
and on the system under study. We will return to this question in Section IIVI 



A. The second coefficient of the fugacity expansion 

The elastic interaction between two particles is decomposed into rotation-invariant sec- 
tors, so that each part depends only on their relative distance r (central interaction): 

V'^ir) =H''- (15) 

^ An iV-particle connected diagram is a graphic representation of N balls linked directly or indirectly by 
lines which represent the correlations coming from interactions or statistical effects. 
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index a standing for the other types of invariance: spin, isospin, charge conjugation, etc.^ 
The scattering matrix 5*1 depends then only on the energy and the angle between the initial 
and final momenta. In this angular momentum representation, the S2 operator can be 
written solely in terms of the phase shifts S^{k): 



{k'l'm'l \klm) 



(16) 



The symmetrization performed by operator A in flTTl) must account for the complete state 
associated with a. It is possible to transfer this symmetrization instruction to the index a. 
Once this is done, the cluster integral of the two-particle system is written as 



62 - i>2 



(0) 



„a 
V 



dE 



'dE 



(17) 



where counts the degeneracy degree of the duly symmetrized states. 

If we carry out a coordinate transformation to the center of mass, the state will behave 
just like a free state (plane wave), that is, it will obey the free Hamiltonian. In order to 
explore this fact, we remember some properties of the relativistic invariant s: 



El + El + - (pi 



P2) 



E' -P' 



(18) 



mi and m2 are the masses of the two particles, E = Ei + E2 is the sum of their energy, and 
the momentum of the center of mass is P = pi + p2- Since E and P are the energy and 
the momentum of the two-particle cluster, u may be understood as the mass of this cluster 
and encapsulates all the information about the interaction. As the interaction calculated 
by ( fTTI) does not depend on the coordinates of the center of mass, the integration and the 
differential operator with respect to E are reexpressed as 

d^P r duj 



E 



92_ 

V 



d^R 



i2n) 



Airi 



S^2 



(19) 



The next step is to integrate on d^P: 



do) 



92 



uj^K2{pu)—Tr 
4m 



M 



1 d - 
5?) —S^ 



du. 



(20) 



For instance, proton-neutron interaction (pn) related to total spin S, total angular momentum J and 
orbital angular momentum L can be represented by a central interaction operator such as VTs+ij^ (r). 
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where M = mi + m2 is the minimum value assumed by u; K2{f3uj) is the modified Bessel 
function. Using f|T6l) . the trace in the angular momentum representation is 



Tr 



-1 d - 

na \ na 
^2 ) 



Equation fl2U|) becomes: 



(21) 



1 f 

h - 4°^ = ^ E E ^2 (2^ + 1) y -^A^2(/5, 



a i=0 



M 



This equation, which reduces to the Beth - Uhlenbeck 



27 



(22) 



28l ] case in the non-relativistic 



limit, is the fundamental tool to obtain the EOS for the pre-nucleosynthesis universe. 



IV. EQUATIONS OF STATE FOR THE PRE-NUCLEOSYNTHESIS UNIVERSE 

The EOS to be built here are supposedly realistic because they account for the interac- 
tions. Nevertheless, we will not actually consider all the four fundamental interactions. 

The gravitational interaction is accounted for only through Einstein (more specifically, 
Friedmann) equations. We will neglect any possible change in the description of statistical 
mechanics due to the curved background of the cosmic manifold.^ 

The electromagnetic interaction, though responsible for the thermalization of the charged 
particles, will not be relevant to the primeval EOS, since: (i) the shielding effect due to the 
existence of opposite charges makes it possible to consider the effective interaction as of short 
range, and the fluid as neutral; (ii) the coupling constant of the electromagnetic interaction 
is relatively small (e^ ~ 1/137). The mean interaction energy between charged particles 
is nearly two orders of magnitude lesser than their mean kinetic energies |l|. The weak 
interaction, though responsible for the thermalization of the neutral particles, is of short- 



range, and its intensity is much smaller [36|] than the strong interaction. It will be neglected. 
Due to its high coupling constant, and despite its short-range, the strong ("hadronic") 
interaction is dominant in the PNS period. As a matter of fact, it is the only one that will 
be taken into account (see Appendix [Bl). 

Justifled by these considerations, we divide the relevant particles for the PNS in three 
categories: 



Ref. 



33| discusses this subject for the dc Sitter solution. 
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1. Ideal ultrarelativistic particles: 7, e~, e"*", fe, z/^, t'e, t'^; 

2. Ideal relativistic particles: fi~,fi^; 

3. Interacting relativistic particles: tt^, tt", 7r°, K~ , K^, A'°, p, p,n,n. 

The equations of state for such a system of particles is the sum of the contributions 
of each species to quantities p and p, keeping in mind what has been said in Section IIIII 
whenever interactions are important. Both p and p can be written as functions of {z,T) or 
in terms of (n, T) - Eqs. fll3m4p . Thus, there are four pairs {p, p} given by the combination 
oi p{z,T), p{n,T), p{z,T) and p{n,T). We have to decide which combination is the most 
suitable. 

There are strong theoretical arguments (based on the series convergences) and outstand- 



ing experimental indications [8|, |37| that the virial expansion p(n, T) is the most convenient 
form for the pressure. The choice between p{z,T) and p{n,T) is more subtle. We shall keep 
the form p{z, T) since it maintains the probabilistic notion inherited by statistical mechan- 
ics. In fact, the energy density is a mean weighted by the Boltzmann factor, and in the case 
of the grand canonical ensemble it reads 

n - — - ^ V ^p-P^-z^- 

p-v-eiv,T,z)^v' ■ ^'^^ 

Comparing this equation to fllUI) and 0131) one sees that only p{z, T) preserves the proba- 
bilistic character order by order. 

Hence, the EOS for the relativistic particles (interacting or not) will be formed by the 
virial series for the pressure and the fugacity series for the energy density: {p{n, T), p{z, T)}. 
Only the first terms of the expansions will be considered, as only two-by-two interactions can 
be calculated or experimentally measured. Following the classification above, the pressure 
and energy density will be 

P = PHad +PR+ Pur = PHad + Pr + (24) 

p = PHad + Pr + Pur, (25) 

where the labels Had, R and UR represents the hadrons, the ideal relativistic particles and 
the ideal ultrarelativistic particles. 



10 



A. EOS for the hadrons 



Pions 7r+, vr", vr'', kaons ^ K~ , K^, and nucleons N = p,p,n,n exhibit spin, isospin 
and charge conjugation symmetries under strong interactions (Appendix [B]). For this reason, 
the set of hadrons can be treated as an interacting system of three components - tt, K and 
N - whose EOS are derived from (]A3cP and (lA4aP . The energy density will be 



PHad{z^,ZK,ZN,kT) = (kT) 



+b2TrnZl + b2KKZK + b2NNZ% + 
+b2TTKZnZK + b2-KNZ-KZN + b2KNZKZN 



(26) 



where the dot indicates differentiation with respect to kT. To this point we have not used 
the fact (already discussed) that /i^ ^ hk — /i at — during the pre-nucleosynthesis period, 
and consequently 2^ 2^ ^i^- ~ 27V — 1- The pressure is 



PHad{nn,nK,nN,kT) = kT {ai^n^ + aixUK + aiNnN+ 

222 

+a27TKnnnK + a27rArn^n7V + a2KN'nKnN) , (27) 



with the numerical densities given by 

nn{Zn,ZK,ZN,kT) = bi^Zn + 2b2jrnZl + b2nKZjrZK + b 2nNZnZN, (28) 
nK{Zn,ZK,ZN,kT) = biKZK + '2b2KKZK + b2nKZ-,rZK + b2KNZKZN, (29) 
nN{Zn,ZK,ZN,kT) = biNZN + 2b2NNZN + b27TNZTrZN + b2KNZKZN- (30) 

The explicit form of the ideal and interaction-related terms are given in Appendix |Al The 
rest mass adopted are: = 0.1396 GeV, rriK = 0.4957 GeV and m^v = 0.93826 GeV . 

In principle, the term 62 should contain all the interaction processes involving the two 
particles to which it refers. That is, it should describe the elastic and inelastic scatterings 
as well as the bound states. Actually, it is possible to argue that, in the construction of the 
EOS for the pre-nucleosynthesis period, bound states and inelastic processes are negligible 
when compared to the elastic channels. 

Taking into account just the elastic processes, the six interaction terms (tttt, vri^, vrA^, 
KK., KN and A^A^) are given by equations of type fl22|) . The explicit expression of the cluster 
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integrals are obtained with the help of all the phase-shift data sets - see Appendix [B] for an 
example. It is necessary to perform an integration involving the modified Bessel function K2 
and derivatives of the phase-shifts. The six integrations have been done numerically with 
the software Mathematica 5.0. Before that, we carefully put back the constants c and % in 
such a way that the length and energy units were fm and GeV , respectively. 

The higher limits of the integrals are different for each 62, since the phase-shift data sets 
are obtained in different energy intervals. Still, most of them lie in the interval from 1 
GeV to 2.5 GeV . As the relevant temperature values for our studies are in the range 20 
MeV < kT < 300 MeV, the chosen values for the higher limits are good enough to include 
all the main contribution from the integral kernels. Besides, the Bessel function K2 assures 
the fast decrease of the kernel values. We consider only angular momenta / = 0, 1 and 2, 
i.e., the S , P and D contributions (Appendix [B|. Once the coefficients bi and 62 have been 
calculated, it is straightforward to obtain the numerical densities (1281) . (!29l) and (130|) . and 
the equations of state for the hadrons, Eqs. (l26l) and ( l27l) . 

The approach described above is not the only one for building hadronic EOS at few 
lundreds of MeV. There is also the treatment called Hadron Resonance Gas (HRG) 



38, m 



40 



4l| describing the Fireballs created in the Heavy Ion Colliders (SPS - Super 



Proton Synchroton - and RHIC - Relativistic Heavy Ion Collider). This technique agrees 
qualitatively with the Lattice QCD calculations [4II and it reproduces the abundance of the 
observed particles relatively well 39| . The HRG modeling is equivalent to ours as long as the 
ideal terms are concerned, but the procedure of including the interaction is rather different. 



The Hadron Resonance Gas accounts for interaction by introducing a excludec 
a la Van der Walls 



42| and by considering the resonance contributions 



38 



volume term 



39|.^ On the 



other hand, our work accounts for interactions through the phase shifts. Because of this, 
other effects - besides the ones related to resonances and hard cores (excluded volumes) - 
are considered. Reference 4J] analyzes these features in detail for the case of A resonance 
in the ttN interaction. It is worth to emphasize that in our approach the resonances are 
only characteristics of the elastic scattering between stable hadrons; the resonances are not 
taken as particles. 



^ The HRG model is similar to the Hagedorn's statistical bootstrap model 43[ in the aspect that both cases 
include the hadronic interactions primary through the resonances. 
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In connection with what we said in the Introduction: the construction of hadronic EOS 
through the S-Matrix was already implemented in Refs. 



29 



30| . However, the only processes 



considered there (e.g. 29|]) are those including the nucleons (scattering NN). The present 
work is more constructive and complete in this sense: all the important statistical quantities 
- pressure, energy density and numerical density - are obtained from the partition function 
which accounts for all the relevant processes of interaction two by two (scattering vrvr, nK, 
vriV, KK, KN and NN). 



B. EOS for the ideal particles 



The ideal (non-interacting) particles are divided into two categories: ultrarelativistic 
(7, e~, e"*", z/g, t'^j, i^e, t'/x) and relativistic (/i~,/i^). The ultrarelativistic sector is relatively 
simple, and has been discussed in Section [III The calculation for the relativistic sector is 
analogous to what we have done for the hadrons. The energy density, for instance, is 



m .2 



The }){^ and 62^^ are the ideal terms determined by ( ]A8a[) and ( lA8bl) . with 
GeV . As in the hadronic case, we set 2;^ — 1 during the PNS. 1 

PR{n^, kT) = kT [ai^n^ + 02^^^^] with n^{z^, kT) 



(31) 
0.10566 

(32) 



C. The complete EOS 

The complete EOS for the energy density p is obtained by substituting Eqs. (|2Bl) . (jl]) and 
(13T]) into (1251) . The plot of p{kT) is shown in Figured] The three curves exhibit similar 
features: all are positive and increase monotonically. The ultrarelativistic components dom- 
inate up to ~ 0.12 GeV. From 0.12 GeV to 0.16 GeV, the other particles (mainly vr and 
p) become important. At 0.275 GeV the ultrarelativistic, ideal-relativistic, and interacting 
sectors correspond respectively to 25%, 33% and 42% of the total energy density. 

The complete EOS for the pressure p comes from the substitution of Eqs.( 1271) . (jlj) and 
fl32|) into flMj) . Figure [2] shows the plot of function p = p{kT). 

It can be seen in Figure [2] that, up to 0.14 GeV , the three curves for the pressure 
present similar behavior: as in the energy density case, all are positive-valued and increase 
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FIG. 1: Plot of the energy density p in GeV/ fm^ as a function of the thermal energy kT given in 
GeV. Three curves are presented: the complete energy density which includes the effects of all the 
relevant particles and their interactions (full line); energy density of all particles (ultrarelativistic, 
relativistic and hadronic) without considering interactions, i.e., all the particles are taken as ideal 
ones (dashed line); energy density for the ultrarelativistic particles (dotted line). 

monotonically. From this point on the interaction contributions grow, and there is a sudden 
change in the slope for the complete system. The curve mitigates its increase rate, reaches 
a maximum at 0.175 GeV, and then starts a steep fall. It becomes negative at around 0.2 
GeV and continues to decrease at a high rate. A comparative graphic analysis of the effects 
on p and p coming from the different kinds of interactions and scatterings reveals that the 
relevant processes for the determination of 62 in scales of kT < 300 MeV are the scatterings 
7171, 7tK and vrA^ related to the S and P partial waves. Therefore, it is fair enough to take 
into account only / = 0, 1 and 2 when calculating 62, just as we did. 

A last point we would like to discuss is the validity range of the expansions p(kT) and 
p{kT). This validity is threatened by three factors: (I) an eventual phase transition due 
to instabilities in the function p{kT); (II) the non-convergence of the pressure and energy 
density series, which would not allow the perturbative approach adopted here; (III) the 
occurrence of the deconfinement of hadrons in quarks and gluons, which would change the 
nature of the particles. Let us comment on these issues: 

(I) There are two necessary and sufficient conditions that assure the stability of the usual 



14 




kJ{GeV) 



FIG. 2: Plot of the pressure p in GeV/ fm^ as a function of the energy kT given in GeV . Interactions 
are accounted for only in the complete case (full line). The conclusion is: interaction drags the 
pressure to negative values, pointing to cosmic acceleration. 



thermodynamical systems 45|], namely: 



T\dT 



> and n~ 



dp 
dui 



> 0. 



(33) 



* / r 



with i = 71, K, N or /i. 

The function p depends only on kT [zi = 1 for every species i), while p depends on kT 
and also on the densities Ui of the four types of particles. Thus, the second condition in fl33|) 
splits into four others, all required to guarantee the system stability. From Figure [H it is 
easily seen that the curve p{kT, {zi}) for the complete set always satisfies fl33|) .^ 

Notwithstanding, the same does not occur with p{kT, {n-i}). Indeed, as the numerical 
densities are positive and increasing functions of kT, Figure [2] makes clear that the pressure 
does not obey condition (1331) : i.e., for some energy value higher than 0.14 GeV the EOS for 
the pressure becomes unstable. In principle, this instability is a strong indication of a phase 
transition and points to a breakdown of the EOS validity. However, this argument is based 
on the usual situation of a thermodynamical system: a gas confined in a recipient of finite 
volume. By its very nature, the universe cannot be assumed to have the same characteristics 
of such a simple and controlled environment. First: the universe has no frontiers. In this 



^ {zi} represents the set of the four fugacities, all of which equals one. 
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case, the surface pressure (which is different from the internal pressure) does not exists. 
And in many cases are precisely the surface effects which generate the phase transitions in 
an ordinary thermodynamical system. Second: in the context of cosmology, the (internal) 
pressure is a direct source of the gravitational field 8|. This particular feature produces 
rather different effects from those engendered by the pressure in usual thermal systems, for 
which one would expect the reduction in the volume as the internal pressure increases. This 
notion, spelled by the second relation in fl33l) . is of great importance, since it prevents the 
ordinary system from disappearing: if dp/dV were positive, nothing would avoid the collapse 
of the gas. On the other hand, in the context of cosmology, an increase of pressure increases 
the gravitational attractive effect, and in this sense, a cosmological system with dp/dV < 
is unstable, i.e., if the universe is not expanding it tends to collapse. We want to emphasize 
the distinction between the effects of the internal pressure in the two types of systems: in 
ordinary equilibrium thermodynamics, an increase of p leads to an increase of the volume V 
of the recipient containing the gas because the shocks of the particles against the walls are 
more and more frequent; but in cosmology the increase of p strengthens the gravitational 
field and produces a net tendency to the reduction of the volume. Hence, it is reasonable 
to say that the stability criterion fl33|l cannot be directly applied in cosmology, and cannot 
be used to rule out the above EOS. From the microscopic point of view, the energy density 
fluctuations 6p grow fast in regions where dp/ dp = < {cs is named the sound velocity in 
the media) and the system becomes unstable. For the usual thermodynamical systems this 
strongly suggests phase transition. Our system, the standard cosmological model, is not of 
this type, though. The universe expands in a rate determined by the relation between p and 
p. This complicates the analyzes and we can not affirm that < necessarily leads to a 
phase transition of the system. The correct treatment of the evolution of 6p should be done 
along with the perturbations in the space-time geometry - FRW metrics. We leave these 
investigations for the future. 

(II) The convergence of the expansions for p and p is related with the formal aspect of 
validity of these series and the correction of the perturbative treatment. This issue is far 
from trivial, since is not possible to sum all the terms for any realistic system in a pure 
thermodynamical context, let alone in the cosmological framework. The virial coefficients 
a, or the cluster integrals b, depend on the sum of the connected diagrams representing 
mutual interactions. And it is extremely difficult to foresee the behavior of the result when 
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particles interact. So, the answer to the question concerning the convergence of p and 
p is inaccessible. Nevertheless, the authors proposed elsewhere js] a toy-model for the 
interactions in the early universe, computed the perturbed EOS (until third order) and 
showed that there is a good indication of convergence for the pressure series. So, we will 
assume in this work the validity of truncating the equations for p, p and n in the second 
order terms. 

(Ill) The QCD coupling constant diminishes as the energy kT increases. So, one could 
expect that at a certain critical temperature kTc, its value becomes sufficiently low to allow 
for the deconfinement of the hadronic matter, giving rise to a system composed by quarks and 
gluons. If such a transition takes place during the PNS, our model ceases to be valid beyond 



kTc'. the EOS have been calculated assuming that the fundament a 
quarks and gluons. According to the lattice QCD calculations 



particles are hadrons, not 



46l |. the critical temperature 



is situated between 150 and 180 MeV. From an experimental perspective, recent results 
from the RHIC 47l| indicate that a sort of transition occurs at these temperature, but it is 
not the deconfinement in a quark-gluon plasma 48|. The experiments run in RHIC found 
what seems to be a new state for the nuclear matter. This state would correspond to a 
perfect fiuid (without viscosity) identified as a CGC {Color Glass Condensate) 49|]. Unlike 
the quark-gluon plasma, the CGC has non negligible correlations, i.e., the nuclear interaction 
processes are relevant. 

There is no doubt left by the experiments that the hadronic matter in the heavy ion 
colliders undergoes a phase transition in energies around 170 MeV. But it is possible that 
this phenomenon does not occur in the primordial universe (or that it occurs in different 
energy scales), due to the various differences between the laboratory and the early universe. 
Specifically, the temporal scale of the events in the possible cosmological QCD transition 
(10~^s) is quite different from the scale of the transition in the accelerators (10~^^s) 50|. 
Other differences - previously cited - are the non-existence of border in the cosmological 
system and the role of pressure as direct source of gravitation. It may happen that the new 
state of nuclear matter (CGC) discovered in RHIC is affected by the absence of borders or 
by the expansion of the universe. Even though we do not have strong arguments against 
the cosmological QCD transition, we will assume the hypothesis that it does not occur in 
the PNS period. 
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V. COSMOLOGICAL CONSEQUENCES 



Let us admit that the proposed EOS, Eqs. (!24|) and (!25|) . are vahd during all the PNS 
period. A first application to cosmology can done through the usual parametrization 



p{kT) =uj{kT)p{kT). 
Substituting flMl) into the second Friedmann equation ([2]), and neglecting A, we find 

^""^ [l + 3c^(A;T)]p(A;T). 



(34) 



(35) 



a J 3 

As the energy density is a positive and increasing function of /cT, the information about the 
acceleration rate depends on the parametric function uj{kT) solely. Isolating this quantity 
in ( IMI) and using the explicit forms of p{kT) and p{kT), the plot of Figure [3] obtains. It 
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FIG. 3: Plot of w as a function of kT. Two curves are shown: the complete to function, including all 
the particles and relevant interactions (full line); ideal to function, which computes all the particles 
but not the mutual interactions (dashed line). 

shows that the complete and the ideal curves for u are practically superimposed until 0.11 
GeV. From this value on, the complete u curve begins to decrease, becoming negative at 
about 0.2 GeV. On the other hand, the ideal u curve remains positive and nearly constant 
up to 0.3 GeV. The effect of the behavior of the complete u function on ( l35l) is to change 
the universe's expansion rate between 0.1 GeV < kT < 0.3 GeV. As we go back in time, 
the universe passes continuously from a decelerated dynamics (pretty close to a radiation 
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dominated era) to an accelerated stage. Comparison of the two lines in Figure [3] shows that 
this acceleration is due to the hadronic interaction processes. 

The next step is to determine the scale factor a{kT) and the function t{kT) associating 
cosmic time and temperature (or thermal energy). We shall use the conservation equation, 

^ = --b + p), (36) 

da a 

and the first Friedmann equation ([2]). We obtain a{kT) by integrating 0361) : 

a{kT) = apnso exp 



dp 
d(kT') 

Kp + p) 



d ikT') 



(37) 



where the label PnsO indicates that the quantity is calculated at the end of the pre- 
nucleosynthesis period which, in energy, is kTpn.so = 20MeV. The time function t{kT) 
is obtained by the combination of ( l36l) and ([2]) (c is re- inserted): 

'^'^ dp 

Equations fl37j) and fl38|) describe the evolution of the pre-nucleosynthesis universe. A 
careful analysis of these equations reveals that both diverge when p = — p and, according to 
the complete EOS (with interactions included), this happens at kTc ~ 0.2767 GeV (critical 
temperature).^ The scale factor always decreases as kT increases. The decreasing rate ^^^^ 
is not so well-behaved: it varies a lot and, as kT approaches kT^, — ^ — oo. This means 
that the scale factor tends to zero when kT —>■ kT^.. In order to study the time function 
around kT^. we have to particularize for the three values of k: 

1. For K = 0, the time interval |tpnso — ic\ ^ oo when kT kT^, because p{kT) and 
■j^^ are always positive. In a plane space-section universe, the time interval that 
would elapse from the initial singularity till tpnso is infinity, and that singularity is 
never attained. 

2. For K = 1, we see from ([2]) that it exists a value of kT close to kT^. where a = 0. 
This value, indicated by kTp, is close to kTc to 39 decimal digits. Defining as 



Not to be confounded with the deconfinement temperature, also usually called "critical temperature" . 
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the present-day value of the scale factor, and as = a{kTB) as the minimum value of 
the scale factor, then it is found that as — 1-7 x lO'^'^ao-^ Thus, in an universe of 
spherical space-section, the initial singularity does not exist; the universe reaches a 



minimum size measured by a^. Recall that recent cosmological data 



5l| slightly favor 



K = 1 (as ^2^0 — 1-003). Combining this type of solution that presents a minimum 



size to the requirement that the closed model presents a.. 



so a maximum size, we obtain 



the so called eternal universe, or bouncing universe 152 1. The time interval from the 
minimum radius till the end of the PNS period is determined by using (138|) and the 
values of Hq and VLto] it is \tpnsQ — — 2.28 ms. 

3. For K = —1, the time interval \tpnso ~ tc\ — ^ 2.28 ms when kT — s> kTc. Therefore, for 
a hyperbolic space-section, the initial singularity is reached in a finite time. 

Notice that despite the different possible evolutions, the universe always presents a max- 
imum temperature kTc. 

The statistical bootstrap model by Hagedorn also predicts a maximum temperature kT^ 
for the universe; but the mechanisms that lead to it are distinct from the causes for our 
critical temperature kTc. In the context of Hagedorn's model, the increase in the total 
energy E is responsible for the rise in the kinetic energy and also for the increase in the 
number of kinds of particles in the system. As E rises it is more advantageous to produce 
new particles than to increase the temperature of the system. This ultimately results in 
an infinite limit for the energy {E oo) and the pressure (p oo) at a finite value of 
the temperature {kTn < oo). On the other hand, kTc is the temperature associated to the 
equality p = — p in our cosmological model with interacting particles - see comments below 
Eq. ( !38l) : kTc appears due to the structure of the equations ([2]) and ( !36l) . 



A. The inflationary regime 

Function uj{kT) evolves continuously from 1/3 (radiation era; kT < 0.02 GeV) and tends 
to —1 (as kT kTc). The inflationary regime (primeval accelerated expansion period) 
occurs when uj{kT) < —1/3. On the other hand, our model requires that uj{kT) > —1. 

^ This calculation considers matter dominance in the interval fcT-yo — 2 x < kT^ < 1 eV and 

radiation dominance in 1 eV < kT^ < 20 MeV , with the initial condition a{kTjo) — oo- 
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We match these requirements by defining the infiationary period as that for which — 1 < 
LoikT) < — |. According to Figure [3|, this corresponds to the energy interval 

0.2356 GeV - kTaccei < kT < kT^ ~ 0.2767GeK (39) 

An early acceleration is required to solve some cosmological problems - horizon, fiatness, 
origin of the inhomogeneities, etc. - without imposing specific initial conditions [s], [lo|, 
An accelerated regime, however, must exhibit some special features to actually rule out 
these problems. We now show that our model indeed eliminates the horizon and flatness 
problems. 

The horizon problem is the lack of causal connection between regions far apart which 
nevertheless exhibit similar physical characteristics. The region of causal connection is 
quantified by the Hubble radius R = {aH)~^, the maximal distance particles can travel 
during an e-fold increase of the scale factor, an increase by the factor e ^ 2.78 4]- The stan- 
dard Bib Bang model tells us that, from the time tpnso until the present day, the maximum 
commoving distance between two causal connected regions grew 9 orders of magnitude: 

Rq 0,PnsoHpnsO ,,,1^9 



~ 1.4 X 10^ (40) 

Let kThor be the energy value above which all the universe's content is in causal contact. In 
order to solve the horizon problem in the context of our model, we need to show that there 
is a kThor < kTc such that 

— = 7.2 X 10 . (41) 

J^hor UPusO-n PnsO 

This is done by using Eqs.([2]) and (1571) . Indeed, for all the values of k, the kThor satisfying 
dS]) respects (kT^ - kThor) ^ 10"^^ 

Analysis of the Hubble radius shows that, from the beginning of the deceleration (kTaccei — 
0.2356 GeV) till nowadays, R has increased 1.75 x 10^*^. But from kThor to kTaccei the Hubble 
radius diminished 1.8 x 10^*^. And for values of time smaller than that corresponding to kThor 
all the universe was in causal contact, allowing the thermalization of the cosmic fiuid. 

The flatness problem may be treated quantitatively through the first Friedmann equation 

According to recent data Qto — 1.003. Now, aH decreases with time if the universe is of 
the radiation-dominated or matter-dominated types (the usual cases). This means that 
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increases with time in radiation and matter-dominated models. If this is so, in the early 
universe ^ lO^^ and we are led to the question: Why was the total density value p 
so close to the critical density pc in the initial instants of the universe when /t = ±1? This 
necessary fine-tuning in the initial condition for p is the fiatness problem. 

The standard model establishes that the energy density at the end of the pre- 
nucleosynthesis period ppnso is close to pc within a precision of 21 decimal digits. Hence, to 
solve the fiatness problem one needs to show that there is a value kTfiat < kTc such that 



la 



1 c y, ]r]~21„2 tt2 



The value kTfiat verifying fH3l) is found to be {kT^ — kTfiat 



10' 



(43) 



10~35 (for K = ±1). From 
kTfiat to the end of the PNS, the scale factor increases 12 orders of magnitude. And this 
large variation happens in a time interval of 2.25 ms, by Eq. fl38|) . The duration of the 
accelerated expansion is much smaller: 0.26 ms. This makes clear that the high infiationary 
rate - a{kT) increase of 10 orders of magnitude in 0.26 ms - rapidly fiattens the universe. 

The values (of the thermal energy, of the scale factor) characterizing the PNS period in 
what concerns the infiationary issues are collected in Table 1. 





At{ms) 


a 

O-PnsO 


R 

Rpn30 




k\T-T,\ (Gel/) 


kTpnsO 





1 


1 


1.5 X 10-21 


~ 0.2567 


kTaccel 


1.99 


5.0 X 10-2 


8.0 X 10-2 


9.55 X 10-24 


~ 0.0411 


kTfior 


2.25 


1.9 X 10-^2 


1.44 X 10^ 


10-3 


10-35 


kTfiat 


2.25 


1.9 X 10-12 


1.44 X 10^ 


10-3 


10-35 



Table 1: Resume of the main results concerning the behavior of our model in the PNS 
period and the accelerated regime. kTpnso = 0.020 GeV is the temperature at the end of the 
PNS; kTacei is the temperature below which there is no acceleration; kT^or is the temperature 
at which the causality problem is solved; and kTfiat is the temperature that rules out the 
fine-tuning in p. Definitions: At = tpnso — t and apnso = 1- 



VI. FINAL REMARKS 



This work presents equations of state (EOS) for the early universe which include the 
interactions among the constituent particles. We argue that the dominant processes af- 
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fecting the pre-nucleosyntesis period (PNS) are those involving the strong interaction, the 
only which has been considered. Total EOS accounting for photons, leptons and hadrons 
have been built, using a phenomeno logical description of the hadronic interaction and as- 
suming thermodynamical equilibrium. Assuming the hypothesis of no deconfinement, some 
interesting points have turned up. 

First, the interactions naturally drove the pressure to negative values. The de Sitter 
model has taught us long ago that an exotic equation of state, p = —p, in which the 
pressure is negative, would give rise to an accelerated expanding universe. Therefore, our 
results could explain the primeval acceleration regime, as an alternative scenario to (scalar 
field) inflation. Indeed, our model is able to connect this initial accelerated stage with a 
decelerated expansion of the type expected for a radiation-dominated universe. In addition, 
the acceleration generated in our interacting cosmic fluid has been shown to be intense 
enough to solve the horizon and flatness problems. 

Another noticeable feature has been obtained: for the models with k = or k = 1 the 
initial singularity is avoided, and a natural explanation for the expansion of the universe 
emerges. In fact, either in an eternal universe (compatible with /t = 1) or in a universe with 
a beginning (corresponding to k = 0), the primeval acceleration produced by the strong 
interaction is capable of engendering an expansion which evolves to a decelerated type and 
continue its dynamics. Nevertheless, let us emphasize that these results have been derived 
for a symmetric universe: the quantity of matter is assumed to equal that of antimatter. 

Such effects are actually more general than the results of the specific model proposed. 
Indeed, in order for then to be valid it is enough to introduce in the FRW cosmology a 
continuous parameter A (related or not to interactions) responsible for the passage from 
p{X) = to the exotic EOS p{X) = — p(A) in such a way that ^ > at all times. 
This type of fitting links smoothly the accelerated and decelerated expansion periods via 
the parametric equation p = up in which u varies between 1/3 and —1. This interval 
automatically eliminates the ghost models(ci; < —1) 53|. 

Besides the mechanisms engendering primeval acceleration, there are other important 
subjects related to the pre-nucleosynthesis period to be discussed. Amongst them, we high- 
light two: the generation of the initial perturbations of the energy density content and the 
matter-antimatter asymmetry. The initial perturbations might be generated in the context 
of our model through the introduction of statistical fluctuation processes 5J]; this possibil- 
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ity shall be investigated in the future. The problem of the matter-antimatter asymmetry 
is more involved: as it is well known, the inflation wipes out all the traces of an eventual 
initial asymmetry and this demands the existence of a mechanism driving a post-inflationary 
matter-antimatter asymmetry. Moreover, a necessary condition to produce an excess of par- 
ticles over anti-particles is the system to be out of the thermodynamical equilibrium 55|. 
This exigence prevents our model to describe this asymmetry once all the numerical densities 
are obtained from the statistical mechanics in thermal equilibrium. This limitation could 
be overcome if we consider out-of-equilibrium phase-transition in a post-inflationary stage. 
A mechanism of this kind was proposed in Ref. 32] and will the investigated in another 
opportunity. 

From a wider perspective, this works calls attention to the fact that particle interactions 
are a direct source of gravitation Q]. When one ignores this truth, simplifications result in 
the cosmological models; but these are perhaps too large to enable a proper description of 
the real universe. The role of the fundamental interactions in cosmology opens new paths 
to the study of the universe evolution, and even if our complete EOS are not valid during 
all the PNS period, the central idea may be applied to model other eras. Up to this point, 
we would not dare to affirm that the inclusion of hadronic interactions is more appropriated 
than the scalar fields to generate the infiation. Many issues solved by the infiationary theory 
have been left untouched here. But the interaction scheme is certainly more suitable from 
the theoretical point of view: they have a clear physical interpretation. 



Acknowledgements 

L. G. M. is grateful to Fundacao de Amparo a Pesquisa do Estado de Sao Paulo (FAPESP) 
and Fundacao de Amparo a Pesquisa do Estado do Rio de Janeiro (FAPERJ), Brazil; R. A. 
and R. R. C. (grant 201375/2007-9) are thank to Conselho Nacional de Pesquisas (CNPq), 
Brazil. R. R. C. and L. G. M. also thank Instituto de Fisica Teorica, Universidade Estadual 
Paulista, Brazil, where this work was initiated. Finally, R. R. C. would like to thank Prof. 
V. P. Frolov and Dr. A. Zelnikov for the kind hospitality extended him at University of 
Alberta. 



24 



APPENDIX A: MULTI-COMPONENT SYSTEMS 



A multi-component system is a set with more than one type of particle or conserved 
quantum number. The grand canonical partition function O for a i?-components system is, 
in analogy to 



oo oo 



e{zr,..,ZB,V,T) =J2-J2 QN,,...,NsiV,T)z^\..z^\ (Al) 

Ni=0 Nb=0 



(See Refs. |56|, l57|.) The grand canonical potential f2 as a function of the fugacities (and 
the temperature) is, then. 



oo oo 



n{z^,...,ZB,T) = ^\neiz,,...,ZB,V,T)= ... bN,,...,N,z^' -Zb" . (A2) 

Ni=0 Nb=0 

where h]si^^,,,^NB the i?-component cluster integrals, with 60,. ..,o = {Ni = N2 = ... = 
Nb = 0). 

The construction of the equations of state - pressure p = p{zi, ..,zb,T), energy density 
p = p{zi, .., zb, T) and numerical densities ni{zi, .., zb, T), 77,2(^1, .., zb, T),..., ub^zi, .., zb, T) 
- is done through the usual mapping from statistical mechanics to thermodynamics, namely 

/ rp\ 00 00 

1-f = ^ = 2^ ••• 2^ bN„...,NBZi '..-Zb", (A3a) 

Ni=0 Nb=0 

(f)C) \ 00 00 

Tt) = E - E NM,...,NbZ^'-4^, (A3b) 

piz„...,ZB,T) . kT^ (§) = ikTf ± ... ± ^^^^-4^ (A3c) 

Symbol {zb} is the set of fugacities, from zi to zb- As in the one-component case, these 
is another form for the EOS: that corresponding to the virial expansion, with the pressure 
and the energy density in terms of the numerical densities. It is: 



"^""7^'^^ =£•••£ «.,.....(T)n^..nL^ (A4a) 



kT 

/i=0 «fl=0 



/ rp\ O" 00 

= E - E ..(rX- (A4b) 



/i=0 /s=0 

The virial coefficients ai^^,,,^ig and ci^,...,Ib are determined by the clusters integrals 6Ari,...,Ar 
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For a 2-component case, the first virial coefficients are 



ao,o = 0, ao,i = ai,o = 1, (A5a) 

_ ~^2,0 _ —^0,2 _ —^1,1 
02,0 - -To ' «0,2 - -Jo ' - I J—y [Abb) 



Co,o — 0, 

C2,0 = —0-2,0, Co,2 



_ bi,o _ ^0,1 

Ci.o — 7 — ) Co,i — 7 — , 
Ol,0 Oo,i 



-ail. 



(A6a) 
(A6b) 



The dot ■ indicates differentiation with respect to kT. Equation (12^ for a 2-component 
system reads: 



02,0 — O2 



00,2 — Oq2 



00 „ 



duj 



2,0 



a 1=0 



2mY 



0,2 



a 1=0 



duj 



(A7a) 
(A7b) 
(A7c) 



1,1 



where mx and my are the masses of the first and second components. The ideal terms are 
determined from the free-system js^ as: 

; (A8a) 
(A8b) 

where gx and gy are the number of internal degrees of freedom of the free-systems; the 
upper (lower) signs refer to the Bose- Einstein (Fermi-Dirac) statistics; and, A(j/3) is the 
relativistic thermal wave-length 

(27r)'j/3 





9x 




9y AO) 


1,0 = 


A^(/3)' 


^0,1 = 




(0) _ 


(±l)^?x 


5(0) _ 


{±l)9Y 


2,0 ~ 


2 A3, (2/3)' 


"0,2 - 





47rm^yi^2(j/3mx,y) 



(A9) 



The generalization for more than two components is straightforward. 
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APPENDIX B: HADRONIC SCATTERING 



Elastic hadronic scattering processes can be described by the phase shifts S in the context 
of the partial-wave formalism. In this approach, one uses the spectroscopic classification 
and separates the "elastic" process involving interactions among pions, kaons and nucleons 
in six types. They are: tttt, ttK, ttN, KK, KN and A^A^. There is a set of phase shifts 5 for 
each of such processes. The 5's are directly or indirectly determined from the experimental 
data. We will discuss in detail how this is done in the case of the pion-pion scattering. 

The "elastic" pion-pion scattering depends on the energy, the total orbital angular mo- 
mentum L and the total isospin I. It is then convenient to introduce the notation 5l i for 
each phase shift. As a scattering of identical particles with integral spin, the associated two- 
pions total state (orbital angular momentum state plus isospin state) must be symmetric. 
We will take / = 0, 1, 2 and restrict the analyses to L = 0, 1, 2 (partial waves S, P and D). 
This restriction is not arbitrary: it is imposed by the experimental data that are available. 
Hence, the relevant phase shifts are those in Table B.l. 



Table B.l: 



5-Wave 


P-wave 


D-wave 


(5o,o with g = 1 
5o,2 with g = 5 


5i 1 with g = 9 


(^2,0 with g = 5 
^2,2 with g = 25 



bosonic 
(2/+1) 



..J.. X iiciDc DiiiltD iwi /wi -scattering. The degeneracy ^icgicc y aL-i^wunto iwi tuc 
symmetry and the projections of orbital angular momentum (2L + 1) and of isospin 



The data for the phase shifts 5o,o and 6i^i are found in Ref. 



59, 



601 ] , which give the results 
for the extensively repeated and measured scattering 7r~p n'-K^n. The data for 5n,2; ^2,0 
and ^2,2 59|, l6ll] are determined through modeling based on the Roy equation 68|, l69|]. The 
S-wave data are shown in Figure HI 

The fits for each one of the phase shifts were done using the software Origin 6.1 via 
the polynomial regression method or the non-linear least squares fitting. The experimental 



The quotation marks were included to indicate that, besides the truly elastic processes (tt^ + p ^ tt^ + p), 
we are considering those with charge conjugation (tt^ + p ^ tt'^ + nj . The isospin symmetries / are also 
taken into account. 
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FIG. 4: Experimental data for the 5-wave tttt scattering and the corresponding fitted curves for 
the functions 5l,i (w). 

uncertainties were considered. The best-fit phase shifts are given by the expressions 

5oo [u] = -1.13 + 4.36CJ - O.llcu^ - 4.52w^ + 3.76cJ^ (Bl) 

and 



2 {to) = -0.03 + 1.23CU - 6.67cu2 + 12.62lj^ - 13.39cj^ + S.STcj^ - 3.11 



UJ 



0.48cj^. (B2) 



Analogously, we can obtain the fitted curves 6lj = (uj) for partial waves P and D. The 
functions 6l,i (uj) - such as fIBip and flB2l) - are differentiated with respect to the energy 
and substituted in equations (lA7atlA7c|) for the cluster integrals b which, in turn, are used 
to obtain the p and p equations of state for the PNS universe. 

The other scatterings {irK, nN, KK, KN and A^A^) are treated in the same way. The 
experimental data used in the description of the pion-kaon scattering are found in Ref. 
and they are good enough only to analyze the values L = 0, 1 - 5'-waves and P- waves. In 
the case of the pion-nucleon scattering, the relevant reference is 



62| 



and we study L = 0, 



1, 2 {S, P and Z^-waves). The kaon-kaon scattering data are (indirectly) obtained from 
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65| : these data refer solely to 



[giI . l66l | and [67| using the separable potential formalism 
the S'-wave KK scattering (we suppose that the processes KK and KK are identical, which 
means that the processes are independent of charge conjugation C . The kaon-nucleon phase 



shifts (for S*, P and D-waves) are in Ref. 63|] - once again we admit independence under C. 



The same Ref. 



63( 1 presents the nucleon-nucleon data for 5, P and D-waves and, in these 
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cases, in addition to the C-independence, it is assumed independence on the isospin /: the 
proton-proton phase shifts are very similar to the neutron-neutron ones. 
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